Theory for Photon- Assisted Macroscopic Quantum Tunneling in a Stack of Intrinsic 
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We propose a theory for photon-assisted macroscopic quantum tunneling (MQT) in a stack of 
capacitively-coupled intrinsic Josephson junctions in which the longitudinal Josephson plasma, i.e., 
longitudinal collective phase oscillation modes, is excited. The scheme of energy-level quantization 
in the collective oscillatory states is clarified in the iV-junction system. When the MQT occurs 
from the single-plasmon states excited by microwave irradiation in the multi-photon process to the 
uniform voltage state, our theory predicts that the escape rate is proportional to iV 2 . This result is 
consistent with the recent observation in Bi-2212 intrinsic Josephson junctions. 
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Since the discovery of macroscopic quantum tunneling 
(MQT) in YBCO grain boundary junctions and Bi- 
2212 intrinsic Josephson junctions (IJJs) 0], a renewed 
interest has been aroused on the MQT in the Joseph- 
son effect from viewpoints of both basic physics and 
quantum-devise applications. In Bi-2212 IJJs the ob- 
served crossover temperature at which MQT appears is 
about 1K0, which is 1-order higher than that in conven- 
tional single-junction systems yl 0, and, furthermore, 
the dissipation effect that prevents a quantum tunnel- 
ing is very weak in spite that the order parameter of the 
high-T c cuprates has the d-wave symmetry with nodes. 
This is because the Josephson plasma frequency in IJJs is 
far higher compared with artificially-made S/I/S Joseph- 
son junctions [5j and the transfer-integral along the in- 
direction vanishes for the nodal quasi-particles. These re- 
markable features in IJJs give great promise for quantum- 
devise applications. In addition, MQT in IJJs is expected 
to provide a new physics originating from the atomic- 
scale multi-junction structure. 

Recently, Jin et al. observed MQT in the switching 
events to the uniform voltage state in Bi-2212 IJJs where 
all the junctions are switched into the voltage state[6|. 
The switching rate observed in the multi-photon process 
is greatly enhanced, depending on the number of stacked 
junctions, N. The observed enhancement in the collec- 
tive switching is proportional to TV 2 , which suggests that 
collective motion of the phase differences in IJJs is re- 
sponsible for the quantum tunneling. In IJJs the longi- 
tudinal Josephson plasma is known to exist as the collec- 
tive motion of the phase differences that propagates in 
the stacking direction of the junctions Q. In this paper 
we formulate a theory for MQT of the collective longi- 
tudinal plasma modes. The scheme of quantum energy 
levels for the collective motion of the phase differences 
in capacitively-coupled N intrinsic junctions is clarified. 
Our theory predicts that the switching rate in the multi- 
photon process is proportional to TV 2 . 

Consider a stack of intrinsic Josephson junctions hav- 
ing tiny in-plane area of W ~ l^m 2 . In the absence of 



an external magnetic field the phase of the order parame- 
ter tpt, £ being the layer index, can be considered uniform 
along the in-plane direction. In this case the coupling be- 
tween junctions in equilibrium state is brought about by 
the long-range Coulomb interaction between charges in- 
duced in the superconducting layers with an atomic-scale 
layer thickness as discussed in 0, that is, the coupling- 
originates from the incomplete charge screening between 
junctions. The charge and phase dynamics in this case is 
well described by the Lagrangian y|, 



-Ej[l - cos(tpe - tpt-\ 
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(1) 



where A® is the scalar potential at the £th superconduct- 
ing layer, -E^,£_i is the electric field inside the insulat- 
ing layer with the dielectric constant e between £th and 
{I— l)th superconducting layers, and s and d are, respec- 
tively, the thicknesses of the superconducting and insu- 
lating layers. The Josephson coupling energy in Eq.(l) is 
defined as Ej — hj c /\e*\, j c being the Josephson critical 
current density. The first term including A® in Eq.(l) 
corresponds to the charging energy which leads to the 
finite charge compressibility To see this we note 

that the canonical momentum which is conjugate to tpg 



is given as 



Pi 



dC h s , . n h . \ 



where ni is understood to be the density fluctuation of 
the Cooper-pairs in the £th layer. Then, the first term in 
Eq.(l) increases the free energy by the amount of oc n 2 , 
which leads to the finite charge compressibility discussed 
by van der Marel and Tsvetkov |9(. The effect of the 
charge compressibility has been observed in optical prop- 
erties of high-T c superconductors |ToL Hl| . 
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The Lagrangian (1) yields the Hamiltonian of the 
form[3, 



Fourier components as 



f E c r 2im i 

iv i c Vri c \l / Q s U= ) \ -pr 1 + 2a(l - cos— — ) \u n u- n 



tm t 

where Qg is the total charge of the iih layer, i.e., Qg = 
e*Wsng, and is the inverse capacitance matrix. Note 
that the condition d 2 <C W is fulfilled even in the intrinsic 
Josephson junctions with a small in-plane area of W ~ 
1/Ltm 2 , since d ~ 1.2nm. From this fact one understands 
that the electric field generated in the junctions is well 
confined inside the stack of the junctions and its direction 
is perpendicular to the junctions, and as a result, the ID 
Coulomb potential is realized almost completely. Then, 
the inverse mutual capacitance diverges linearly as 9] 



And 

eW ' 



for <>m. 



(4) 



This result indicates that the long-range nature of the ca- 
pacitive coupling between junctions should be correctly 
incorporated, that is, any truncation of the inverse capac- 
itance matrix in Eq. (3) cannot be accepted in the intrin- 
sic Josephson junctions. As shown in [12j, the Hamilto- 
nian (3) can be transformed into more tractable form by 
the canonical transformation. We choose (6g,ug) defined 
as 



h = ¥>e - ipe-i, 



Ug 



(5) 



as the canonical variables instead of (ipg,pg). It is easy to 
see that the canonical commutation relation [6g,u m ] — 
ifi8i m holds if [tpg,p m ] = ihSg m . The Hamiltonian (3) 
is greatly simplified if one uses the canonical variables 
(B e ,u t ) as 
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n = XXfr [(! + 2a ) u ? ~ 2cm ^+i] + E J [! - cos h] 

1=1 

(6) 

with E c = 2nde* 2 /We and a = efi 2 /sd. Note that the 
Josephson coupling term is diagonal and the coupling be- 
tween junctions survives only for nearest neighbors. The 
coupling constant a is identical with the one introduced 
in the capacitively-coupled classical intrinsic Josephson 
j unctions y], la] ■ In Bi-2212 a takes a value of ~ 0.1 — 0.2, 
depending on the doping level [§1 Il3l IT^|. Let us discuss 
the small quantum phase oscillations in the iV-junction 
system. We consider the case of N ^> 1 and neglect 
the boundary effect in the following calculations. For 
small phase oscillations one can use the approximation, 
1 — co& 6 1 ~ \®\- Using the periodic boundary con- 
dition, we express ug and Qg in terms of the Fourier 
series expansions as ug = n ti„ exp(i2irn£ / N) , 

9 e = AT -1 / 2 ^ n n exp(i27rn^/JV'), where the Fourier 
components (u n -,9n) satisfy the commutation relation, 
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The Hamilotonian (6) in the har- 
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(7) 



where n c is given as N — 2n c + 1. Then, in terms of the 
creation and annihilation operators c n and c n defined as 
n = (c„ + cL„)/(\/2if„) and u n = -iHK n (c n - ct n )/\/2 
with K n = (Ej/ujn) 1 / 2 , one can diagonalize the Hamito- 
nian as 



n = n. 



n=—n 



< t l s 
Un(C n C n + - ), 



where 



2nn s 

l + 2a(l-cos— ), 



(8) 



(9) 



monic approximation can be expressed in terms of the 



with Lu p i being the Josephson-plasma frequency defined 
as hujpi = y/2E c Ej. This phase oscillation mode is identi- 
fied with the quantum version of the longitudinal Joseph- 
son plasma. The dispersion relation given in Eq.(9) is the 
same as the classical one 0,0- The longitudinal Joseph- 
son plasma has been observed in the microwave absorp- 
tion experiments in Bi-2212^^ Its frequency u) p \ 
is located in the range of a few hundreds GHz at T = 
OK^I- As seen from Eq.(9), we have N plasma modes in 
the iV-junction system, which is crucially different from 
the single-j unction system, and their frequencies range 
from oj n= a = oj p \ up to w„ =±A r/2 = \/l + 4aw p i. In the 
case of a ~ 0.1 (under-doped Bi-2212 case) the relation, 
VI + 4ao-> p i — uj p \ <C ojp\, holds. Then, the lowest single- 
plasmon state is well above the highest zero-plasmon 
state (see Fig.l). From this observation one expects 
the energy scheme as depicted schematically in Fig.l 
for the capacitively-coupled Bi-2212 intrinsic Josephson 
junctions in the quantum regime, if the harmonic approx- 
imation for the phase oscillation modes is assumed to be 
valid up to the single-plasmon state. In this figure we 
also describe the washboard-like effective potential which 
confines the collective phase oscillation modes. The tilt 
of the potential can be generated by a bias current. Note 
that in the presence of a bias current / the plasma fre- 
quency is shifted as u> p i — > u> p i(I) = uj p \(1 — I 2 ) 1 / 4 with 
/ = I/j c by the effect of the coupling term —EjIJ^g ®t- 
If /pi = ^ P i/27T ~ 100GHz, we have the level spacings 
hio p i ~ 4.8K between zero-plasmon and single-plasmon 
states and hu} n+ i—hio n ~ 0.23K between adjacent plasma 
modes for a — 0.1 and N = 50. From this estimation the 
quantum effect originating from the discreteness of the 
n-plasmon states may be expected at T ~ IK. The level 
splitting between the plasma modes can be seen at tem- 
peratures below ~ lOOmK. 




FIG. 1: Energy-level quantization of the longitudinal collec- 
tive phase oscillation modes in the presence of a bias current. 
In the iV-junction system single- plasmon states are composed 
of N oscillation modes which are nearly degenerate in the case 
of a <C 1. The zero-plasmon state is formed by the zero-point 
motion of the N plasma modes. 
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FIG. 2: Two-photon process that causes the transition |0 >^ 
\u)k > in the third-order perturbation. In this process the vir- 
tual plasmon of mode k — with energy 2u, which is excited 
in the two-photon process, excites a plasmon of mode k by 
the interaction V3 when 2v = Wfe. In this figure Oo > and 
I Ojfe > denote, respectively, the vacuum states of Oth and fcth 
plasma oscillations, and \ujq > and \u>k > the single- plasmon 
states of these oscillations. 



Let us now study the MQT in the intrinsic Joscphson 
junctions in which the longitudinal Josephson plasma is 
excited. Suppose that the single-plasmon states show the 
MQT to the voltage state. Since the longitudinal plasma 
oscillations are the coherent motion of all the junctions, 
it is reasonable to assume that the voltage state switched 
from the single-plasmon state is uniform, that is, all the 
junctions are in the voltage state. This switching phe- 
nomenon should be discriminated from the MQT to the 
first resistive branch^- The quantum switching to the 
first resistive branch is induced by the MQT of the non- 
linear localized mode such as the discrete breather [18| to 
the localized rotating modedE^. The MQT of the dis- 
crete breather will be discussed in a forthcoming paper. 

In this paper we study the photon-assisted collective 
MQT in the intrinsic Josephson junctions, because the 
MQT to the uniform resistive state has been observed 
only in the systems under the microwave irradiation up 
to now 6] . In the following we focus on the origin of 
the intriguing ./V 2 -dependence of the MQT rate recently 
observed in Bi-2212. For this we consider the following 
third-order anharmonic interaction term that appears in 
the presence of a bias current, 

v 3 = -1ejiJ2o! 

t 

= ~ V/N ^ ^rn^njns, (10) 

V ni+n 2 +"3=0,±JV 

with 7 being some constant. Note that this interac- 
tion term induces the coupling among the plasma modes, 
i.e., the mode-mode coupling. The conservation law in 
Eq.(10), ni + n 2 + n 3 = 0, ±iV, comes from the peri- 
odic boundary condition. In the real systems composed 



of N ~ 10 — 100 junctions this conservation law is un- 
derstood to hold approximately. One expects that this 
anharmonic interaction term becomes important for the 
excited states, i.e., the single-plasmon states, as the bias 
current is increased. 

Let us investigate the excitations of the plasma mode 
by microwave irradiation, that is, the transition from the 
ground state |0 > to fcth single-plasmon state \ujk > (see 
Fig. 2). We assume that the coupling between microwave 
and the phase differences is generated by the oscillating 
current induced by the microwave. Then, the interaction 
with the microwave of frequency v is described by the 
term as 

V M {t) = EjI M sin vt ^2 8e = VNgu(co + cj) sin vt, 
t 

( n ) 

where gu — EjIm/(2V2K ), Im being the intensity of 
the normalized oscillating current. From Eq.(ll) one un- 
derstands that the microwave can excite only the uniform 
plasma mode of n — 0, Iujq >> in the one-photon process. 
The transition rate l/r' 1 ' in this process is easily derived 
as 

^jy - NnglSMi) - v). (12) 

On the other hand, in the multi-photon processes in 
which the mode-mode coupling is incorporated, the 
plasma modes with k =/= can be excited as shown in 
the following. Note that Eq.(10) includes the interaction 
term, 

^3 — 7^7 C l C mfc C L fc : (13) 
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where //. = "/EjI(SKkK^ rik )^ 1 / 2 and mk is the index 
satisfying the relation k + 2m k = 0, ±TV. Let us study 
the third-order perturbative processes caused by the in- 



teraction Hint — Vu{t) + V% + V 3 . The two-photon 
process is included in the third-order perturbation and 
the transition, |0 >^ \uj k >, is possible in the presence 
of V 3 k as seen in Fig.2. Note that V M {T) cx VN and 
V 3 k oc (a/TV) -1 - Hence, one understands that the tran- 
sition matrix element < f\i >=< uj k \U^ (t, — oo)|0 > 
in the two-photon process is proportional to V~N, where 
U^ 3 '(t, — oo) is the third-order time-evolution operator 
given by 

t/ (3) (t,-oo) = (- l -f f atx f 1 dt? ! <!/ 



2 / U(. 3 

oo 



e e{tl+t2+t3) H in t(h)H in t{t2)H in t(t 3 ) , (14) 



and thus the transition rate, r fc = 4r\ < 
Wfe|f/^(t, — oo)|0 > | 2 , is proportional to TV. In fact, 
from the simple calculations we obtain 



i =N _ <9 2 M h) 2 



S(u k -2v). (15) 



Then, if the MQT rate of the single-plasmon state \oj k > 
to the voltage state is given by r k , the photon-assisted 
switching rate is obtained as r k ~ 1 T k . In the MQT at tem- 
peratures around ~ IK the splitting between the plasma 
modes is smeared out for the systems with TV ~ 50. 
Then, we find the total switching rate in the two-photon 
process for v ~ to p i/2 as 



E 

k——n. 



-T k cx TV 



-r , 



(16) 



which is proportional to TV 2 . The TV 2 -dependence of 
the switching rate is also obtained in the three-photon 
process. Note that the transition, |0 >— * \u>k >, in 
the three-photon process is possible in the fifth-order 
perturbative process with respect to 7ii n t- Then, from 
Eqs.(10) and (11) we expect the relation, < f\i >cx 
(\/TV) 3 x (\/N)~ 2 = v^V, for the three-photon process 



in the fifth-order perturbation, which is the same as in 
the two-photon process. This observation indicates that 
the total switching rate is also proportional to N 2 in the 
three-photon process. Thus, one may conclude that the 
MQT rate to the uniform resistive state is proportional to 
TV 2 , if the MQT happens collectively in the multi- photon 
process. This result is consistent with the recent experi- 
mental result in Bi-2212 intrinsic Josephosn junctions^]. 
Note that the A f2 -dependence dose not appear in the 
switching rate by way of the single-photon process. This 
is because the single-plasmon state of fc ^ cannot be 
available in the single-photon process. Thus, our theory 
predicts that the switching rate in the single-photon pro- 
cess is proportional to TV. Furthermore, from the result 
given in Eqs.(15) and (16) it is also predicted that the 
frequency spectrum of the switching rate, i.e., T^ 2 \v), 
splits into N/2 peaks having centers at v = Lu k /n at low 
temperatures, say T ~ lOmK, in the n-photon process 
if the damping originating from the coupling with the 
quasiparticles or environment is weak enough. 

In summary we have formulated MQT for the collective 
longitudinal Josephson plasma modes in intrinsic Joseph- 
son junctions, assuming that the switching events to 
the uniform voltage state occur from the single-plasmon 
states. The scheme of energy-level quantization for the 
collective oscillation modes has been clarified in the case 
of Bi-2212 IJJs. In the TV-junction system the single- 
plasmon states are composed of N different oscillation 
modes whose frequencies range from u> p i to \/l + 4aw p i. 
These modes can be excited by the microwave irradiation 
in the multi-photon process. Then, there are TV channels 
in the MQT of the collective phase oscillation modes in 
the multi-photon process. We have also shown that the 
transition rate to each single-plasmon state in terms of 
microwave irradiation is proportional to TV. Therefore, 
the photon-assisted switching rate to the uniform voltage 
state is proportional to TV 2 in the multi-photon process 
around IK. 
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